P BEST AVAILABLE COPY 



FRANQAISE 



f 




PCT/FR/20 0 4 / 01 7 ? 



INSTITUT 
NATIONAL DE 
LA PROPRIETE 
1NDUSTRIELLE 




BREVET D. ■ INVENTION 



CERTIFICAT D'UTILITE - CERTIFICAT D'ADDITION 



COP1E OFRCIELLE 

Le Directeur general de I'lnstitut national de la propriete 
industrielle certifie que le document ci-annexe est la copie 
certifiee conforme d'une demande de titre de propriete 
Industrielle deposee a I'lnstitut. 



Fait a Paris, le . 



JHIIff 2004 



Pour le Directeur general de I'lnstitut 
national de la propriete industrielle 
Le Chef du Departement des brevets 



DOCUMENT DE PRIORITE 

PRESENTE OU TRANSMIS 
CONFORMEMENT A LA 
REGLE17.1.a)OUb) 



Martine PLANCHE 



SIEGE 



INSTITUT 26 bis - ™ e de Saint-Petersbourg 

N A T i n ri a i n c 75800 PAR,S ced <* 08 

N A T 1 O N A L D E . 3g ^ 53 04 ^ ^ 

PROPRIETE Telecople:33(0)153 04 45 23 

INDUSTRIELLE www.inpl.fr 




26 bis, rue de Saint Petersburg 
75800 Paris Cedex 08 

Telephone : 01 $3 04 53 04 Telecopie : 01 42 94- 86 54 



REMISE PES PIECES 
DATE 

ueo 11 JUIL 2003 
75 INPI PARIS 

N° 0* ENREQilSTREM ENT 

NATIONAL ATTRIBU* PAR IWI 0308S 1 9 

DATEPH D^POT ATTRIBUTE 
PAR L'INPl 



f t ML 2003 



Vos references pour ce dossier 

(facullalif) 105136/SH/SSD/TPM 

Confirmation d'un ttepdt par tglggopte 
NATURE DE LA DEMAND E 



BREVET D'INVENTION 
CERTIFICAT D'UTILITE 

Code do la proprlete intellectuello - Uvre VI 

REQUETE EN D&JVRANCE 1/2 

Cetimprime est a rempllr Hsiblement a I'encre noire 



N» 11354*01 



[NOM ET ADRESSE DU DEMANDEUR OU DU MANDATAIRE 
A QUI LA CORRESPONDANCE DOIT ETRE ADRESSEE 

" COMPAGNIE FINANCIERE ALCATEL 
Departement PI 
Stephane HEDARCHET 
5, aie Noel Pons 
92734 Nanterre Cedex 



Be5«aw/?6C89» 



Demande de brevet 



□ N a attribue par HNPj a la telecopje 
Cochez I'une des 4 cases sulvantes 



Demande de certificat d 'utllite 
Demande divisionnalre 

Diimdtulede brewl miliak 
ou (lemattdode certified fuimini liak 
Transformation d'une demande de 
brevet europeen Demuiuk de bre&u initiate 
Q TITRE DE L'INVENTION (200 can**™ o^pac.s m ^ mum) 

COMPRESSION CONTEXTUELLE D'IMAGES NUMERIQUES 



m 



□ 

tr 




Date 



DECLARATION DE PRIORITE 
OU REQUETE DU BENEFICE DE 
LA DATE DE DgPOT D'UNE 
DEMANDE ANT^RIEURE FRANQAISE 



DEMANDEUR 



Norn ou denomination sociale 



Pays ou organisation 
Date I / / 



Pays ou organisation 
Date I L L 



Pays ou organisation 
Date I / / 



□ yil y a d'autres priority cochez la case et utilise* fimprime ccSuite* 

□ 8H y a d'autres dematideurs, cochgz la case et utilise* I'ImprTm* «Sufte» 



ALCATEL 





BREVET D'INVENTION 
CERTIFICAT D'UTILITE 

REQUETE EN DEUVRANCE 2/2 



REMISE PES PlfcCES 

^ 11 JUIL2O03 
75 INPl PARIS 

I N° D'ENREGISTREMENT Q3085 1 9 
I NATIONAL ATTRlSUi PAR L'MP> 


0B 540W/26G899 


I Vos references pour ce dossier ; 1 

\(faadlaUf) | 


105136/SH/SSD/TPM ^ 


IQ mandataire 1 




| Norn 


HEDARCHET 


| Pr6noro 


Steohane 


| Cabinet ou Societe I 


Compagnie Financiere Alcatel 


1 N °de pouvoir permanent et/ou 1 
1 de lien contractual J 




1 Adresse 


Rue 


5, rue Noel Pons 


Code postal etville 


92734 1 NANTERRE Cedex 


I N° de telephone (faadlalif) 




1 N° de telecopie (facidtatif) 




1 Adresse electronlque (faadtafij) 




j H 1NVENTEUR (S) 




| Les inventeurs sont les demandeurs 


□Out * 

1 gj Non Dans ce cas fourn ' r une designation d'inventeur(s> separ£e | 


j |3 RAPPORT DE RECHERCHE 


| Uniqueitient pour une demande de brevet (y compris division et transformation) 


1 Etablissement imm§diat 
1 ou etablissement differe 


ho 
□ 


I Paiement echelonnS de la redevance 


I Paiement en trois versements, unlquement pour les personnes physiques 
□Qui 


| M REDUCTION DU TAUX 
1 DES REDEVANCES 


1 Uniquemant pour les personnes physiques 

1 ORequise pour la premfere fois pour cette invention QatndmunavUdenan-imtmilfon) 
1 IZlRequise anterieurement h ce d6pot (jahtdw une copk de In decision d'udmismn 
1 pour cette invention ou mdiquersa rtfSretice) : 




1 Si vous avez utilise Timprimd «Suite», 
1 indiquez le n ombre de pages jointes 


I 




US SIGNATURE dOOOfiXMttl&C 

»M MANDATAIRE Stephane HEDARCHET / LC 40 B 

I ^ 


VISA DE LA PREFECTURE 
OU DE L'lNPI 

L MARIELLO 



La loi n°78-17 du 6 Janvier 1978 relative a I'informatique, aux fichiers et aux liberies s'applique aux rdponses feites k ce formulaire. 
Elte garanttt un droit d'accfcs et de rectification pour les donnSes vous concernant aupriss de I'INPL 



COMPRESSION CONTEXTUELLE D'IMAGES NUMERIQUES 



Invention oonceme le domaine de la compression/decompression 
de donnees d'lmages numeriques par transformee en ondelettes. 

Le codage par transformee en ondelettes consiste classiquement a 
coder une image en trois etapes. Une premiere etape est destinee a 
decorreler I'image en transformant ses pixels en coefficients d'ondelettes 
repartis dans des sous-bandes se situant a plusieurs niveaux de resolution 
Une deuxieme etape est destinee a quantifier les coefficients d'ondelettes 
obtenus lors de la premiere etape. Une troisieme etape est destinee a 
effectuer un codage entropique des coefficients d'ondelettes quantifies lors de 
la deuxieme etape. 

Comme le sait I'homme de I'art, ce type de codage ne peut- etre 
optimal qu'a condition que trois hypotheses soient verifiees : la transformee 
en ondelettes doit etre orthogonale, les coefficients d'ondelettes dolvent etre 
codes par un codage entropique d'ordre 0 dans lequel on n'exploite aucune 
dependance statistique entre coefficients, et on doit se trouver dans les limites 
de la « quantification a haute resolution », dans laquelle les densitel de 
probability des coefficients d'ondelettes varient peu a I'interieur d'une boTte de 
quantification. 

Lorsque ces hypotheses ne sent pas verifiees, il faut en theorie faire 
une allocation de debit complete, en choisissant un pas de quantification a 
pr.ori different pour chaque sous-bande. Dans la pratique, le codage demeure 
qualitativement valable lorsque les hypotheses ne sont plus formeilement 
verrfiees. C'est notamment le cas lorsque la transformee en ondelettes n'est 
pas orthogonale, mais bi-orthogonale, ou lorsque les coefficients d'ondelettes 
d'une meme sous-bande ne sont pas codes independamment du fait qu'ils ne 
sont pas independants. De meme, lorsque I'hypothese de quantification a 
haute resolution n'est pas verifiee, on peut modifier legerement le codeur en 
remplacant son quantificateur uniforme par un quantificateur a intervalle 
central (intervalle de quantification centre en zero ou Zone Morte) de faille 
double. Cela resulte du fait que la distribution des coefficients d'ondelettes 
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dans chaque sous-bande est laplacienne (ou gaussienne generalis6e) avec 
un pic en 0 oO eiie ne peui plus etre raisonnablement consideree comme 
quasi constante surtout aux debits faibles (c'est-&-dire aux forts taux de 
compression). II est ici rappel6 qu'une fonction laplacienne est un cas 
5 particulier d'une gaussienne generalisee. 

Pour tenter d'ameliorer le codage dans les conditions precitees 
(transformation bi-orthogonale ou transformation ne conduisant pas & 
Tindependance de ses coefficients) il a 6te propose de tirer partie de la 
dependance entre coefficients d'ondelettes voisins dans une meme sous- 
10 bande en realisant un codage entropique contextuel prenant en compte une 
densite de probabilite conditionnelle en lieu et place d'une simple densite de 
probabilite. Cette densite conditionnelle modelise la loi d'apparition d'un 
coefficient d'ondelettes conditionnellement a son voisinage. 

Dans ce qui suit, on entendra par « dependance intra-bande » une 
is dependance entre un coefficient q.k.k 1 (ou Wj^[n,m]) et des coefficients voisins 
de meme echelle ou sous-bande Wjk, comme par exemple q,k-i,k' et q^k'-i ou 
Cj,k-i,k' et q,k-i t k'-i- De telles dependances intra-bande sont notamment 
exploitees par les codeurs qui effectuent un codage des positions des 
coefficients significatifs en mode « run-length » ainsi que par les codeurs de 
20 type EBCOT/JPEG2000. 

Par ailleurs, dans ce qui suit, on entendra par «d§pendance inter- 
bandes » une dependance entre un coefficient q^x et un coefficient plac£ a la 
meme position a T^chelle superieure, comme par exemple Cj+i t k/2 ( kV2- De telles 
dependances inter-bandes sont notamment exploitees par les codeurs dits 
25 «en arbre de zeros » (developpes par Shapiro) et le codeur SPIHT 
(developpe par Said et Peariman). 

Ces dependances decoulent de la representation de chaque 
coefficient d'ondelettes quantifie Q(Cj t k,k ,n ) sous la forme d'une chaihe 
d'informations comportant son bit de coefficient significant (N jt k,k ,n = 0 ou 1 ), et 
30 lorsque ce bit vaut un (1), son signe (Sj,k,k ,n qui vaut ±1) et sa valeur absolue 
(Aj, k /). 

Par exemple, le codage par arbre de zeros exploite la dependance 
statistique entre les bits de coefficient significatif (Nj,k,k' n )> tandis que le codage 
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de type EBCOT exploite la dependance statistique entre les bits de coefficient 
significatif (Nye*") et la dependance entre signes S J>k ,k ,n de coefficients 
d'ondelettes voisins quantifies avec differents pas. 

Ces codages contextuels mettent en ceuvre des m£thodes de 
traitement particulierement complexes qui imposent aux codeurs d'apprendre 
autant d'histogrammes qu'il y a de contextes, ce qui limite notablement leur 
Vitesse d'adaptation et interdit leur utilisation dans le cadre de la compression 
a haut d£bit d'entree (images de grande taille) et en temps reel, notamment 
lors d'observations spatiales. 

Uinvention a done pour but d'am61iorer la situation. 
Elle propose a cet effet un precede dedie & la compression de 
donnees damages numeriques par transformee en ondelettes et comprenant 
une §tape consistant a transformer des donnees d'image numeriques en 
coefficients d'ondelettes r6partis dans des sous-bandes, une etape consistant 
& quantifier les coefficients d'ondelettes* et une etape consistant a coder de 
fagon entropique les coefficients d'ondelettes prec§demment quantifies. 

Ce proced§ se caracterise par le fait, d'une part, qu'H comprend, entre 
les 6tapes de transformation et de quantification, une 6tape consistant a 
estimer, pour chaque coefficient de chaque sous-bande d'une image, de : 
premier et second jeux de parametres de prediction, associes respectivement- 
a des directions dites « nord-sud » et « ouest-est », en fonction des valeurs 
des coefficients d'ondelettes de ses voisins nord ou ouest, et d'autre part, 
qu'a Petape de codage entropique, d'une part, on determine pour chaque 
coefficient de chaque sous-bande de I'image des valeurs de prediction de 
l'esp<§rance et de la largeur d'une fonction laplacienne (ou gaussienne 
g§neralisee), representative de sa densite de probability en fonction des 
premier et/ou second jeux de parametres de prediction et du coefficient 
d'ondelettes quantify de son voisin nord ou de son voisin ouest, et d'autre 
part, on code de fagon entropique les coefficients d'ondelettes quantifies a 
I'aide des esperances et largeurs associees determines. 

Grace aux estimations des parametres de prediction, la phase de 
prediction de l'§tape de codage entropique est notablement simplifiee, ce qui 
permet de compresser (ou comprimer) des donnees d'image numeriques a 
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haut debit et en temps reel, notamment lors d'observations spatiales. 

Le procede selon I'invention peut comporter d'autres caracteristiques 
qui pourront Stre prises separement ou en combinaison, et notamment : 

- lorsque cela s'avere possible, on peut determiner lors de I'etape de codage 
s entropique, pour la fonction laplacienne de chaque coefficient d'ondelettes 

quantifie, des valeurs de prediction d'esperance nord et ouest et des 
valeurs de prediction de largeur nord et ouest. La vaieur de prediction 
d'esperance nord (ou ouest) est alors preferentiellement definie par le 
produit entre un premier parametre de prediction du premier (ou second) 

10 jeu et la vaieur du coefficient d'ondelettes quantifie voisin nord (ou ouest), 
et la vaieur de prediction de largeur nord (ou ouest) est definie par la 
somme d'un deuxieme parametre de prediction du premier (ou second) jeu 
et du produit d'un troisieme parametre de prediction du premier (ou 
second) jeu et de la vaieur absolue du coefficient d'ondelettes quantifie 

is voisin nord (ou ouest). Dans ce cas, quatre cas doivent etre envisages 
selon I'environnement du coefficient d'ondelettes concerne. En I'absence 
de coefficient voisin nord, on identifie les valeurs de prediction d'esperance 
et de largeur de la fonction laplacienne aux valeurs de prediction 
d'esperance ouest et de largeur ouest. En I'absence de coefficient voisin 

20 ouest, on identifie les valeurs de prediction d'esperance et de largeur de la 
fonction laplacienne aux valeurs de prediction d'esperance nord et de 
largeur nord. En I'absence de coefficients voisins ouest et nord, on effectue 
un codage laplacien fixe. Enfin, en presence de coefficients voisins ouest 
et nord, on identifie les valeurs de prediction d'esperance et de largeur de 

25 la fonction laplacienne respectivement a la demi-somme des valeurs de 
prediction d'esperance nord et ouest et a la demi-somme des valeurs de 
prediction de largeur nord et ouest, 

- on peut estimer les parametres de prediction des premier et second jeux 
par regression. 

so • Par exemple, on peut estimer les premiers parametres de prediction 
des premier et second jeux par une premiere regression obtenue en 
resolvant une equation de type W[n+1 ,m] = e.W[n,m], ou W[n,m] est la 
vaieur d'un coefficient d'ondelettes et W[n+1,m] est la vaieur du 



coefficient d'ondelettes voisin de W[n,m], et 9 est un premier 
parametre de prediction. On peut par exempie resoudre la premiere 
equation a I'aide de Ja technique des moindres carres, le premier 
parametre de prediction d'un premier ou second jeu etant alors egal au 
rapport entre les covariance et variance non normalisees des valeurs 
des coefficients d'ondelettes associes au premier ou second jeu. 
• De meme, on peut estimer, d'une part, les deuxieme et troisieme 
parametres de prediction du premier jeu par une seconde regression 
obtenue en resolvant un premier systeme d'equations de type 
W[n,m] - 9 N .W[n,m-1] = a N + p N W[n,m-1] , ou W[n,mj est la 
valeur d'un coefficient d'ondelettes, est la valeur du coefficient 

d'ondelettes voisin nord de W n , m , et 6 N , a N et p N sont respectivement 
des premier, deuxieme et troisieme parametres de prediction du 
premier jeu, et d'autre part, les deuxieme et troisieme parametres de 
prediction du second jeu par une seconde regression obtenue en 
resolvant un second systeme d'equations de type W[n,m] - e 0 .W[n- 
1,m] = cc 0 + Po W[n-1,m] , ou W n . 1>m est la valeur du coefficient 
d'ondelettes voisin ouest de W n>m , et Go, ct 0 et p 0 sont respectivement 
des premier, deuxieme et troisieme parametres de prediction du 
second jeu. On peut par exempie resoudre les premier et second 
systemes d'equation par la technique des moindres carres, 

- lorsqu'un coefficient d'ondelettes quantifie voisin presente une valeur nulle, 
on peut le substituer par le produit du pas de quantification, utilise lors de 
I'etape de quantification, et d'une constante choisie, strictement superieure 
a zero et inferieure ou egale a un (1), et plus preferentiellement egale a 1/3, 

- lors de I'etape de codage entropique, on peut proceder a la discretisation 
de la fonction laplacienne apres avoir determine les valeurs de prediction 
de I'esperance et de la largeur. En variante, on peut proceder au codage 
entropique par iteration apres avoir determine les valeurs de prediction de 
I'esperance et de la largeur. 

L'invention porte egalement sur un dispositif dedie a la compression 
de donnees d'images numeriques par transformee en ondelettes et 
comprenant des moyens de traitement munis, d'une premiere part, d'un 
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module charge de transformer les donnees d'image numeriques qu'il recoit en 
coefficients d'ondelettes repartis dans des sous-bandes, d'une deuxieme part, 
d'un module charge de quantifier les coefficients d'ondelettes, et d'une 
troisieme part, d'un module charge de coder de facon entropique les 
5 coefficients d'ondelettes quantifies. 

Ce dispositif se caracterise par le fait, d'une part, que ses moyens de 
traitement comprennent egalement un module charge d'estimer, pour chaque 
coefficient de chaque sous-bande de I'image, de premier et second jeux de 
parametres de prediction, associes respectivement a des directions dites 
10 «nord-sud» et «ouest-est», en fonction des valeurs des coefficients 
d'ondelettes de ses voisins nord ou ouest, et d'autre part, que son module de 
codage entropique est charge, d'une part, de determiner, pour chaque 
coefficient de chaque sous-bande de I'image, des valeurs de prediction de 
I'esperance et de la largeur d'une fonction laplacienne (ou gaussienne 
is generalisee), representative de sa densite de probabilite, en fonction des 
premier etfou second jeux de parametres de prediction et du coefficient 
d'ondelettes quantifie de son voisin nord ou de son voisin ouest, et d'autre 
part, de coder de facon entropique les coefficients d'ondelettes quantifies a 
I'aide des esperances et largeurs associees determinees. 
20 Le dispositif selon I'invention peut comporter d'autres caracteristiques 

qui pourront etre prises separement ou en combinaison, et notamment : 
- son module de codage entropique peut etre agence de maniere a 
determiner pour la fonction laplacienne de chaque coefficient d'ondelettes 
quantifie, lorsque cela s'avere possible, des valeurs de prediction 
25 d'esperance nord et ouest et des valeurs de prediction de largeur nord et 
ouest. La valeur de prediction d'esperance nord (ou ouest) est alors 
preferentiellement d§finie par le produit entre un premier parametre de 
prediction du premier (ou second) jeu et la valeur du coefficient 
d'ondelettes quantifie voisin nord (ou ouest), et la valeur de prediction de 
30 largeur nord (ou ouest) est definie par la somme d'un deuxieme parametre 
de prediction du premier (ou second) jeu et du produit d'un troisieme 
parametre de prediction du premier (ou second) jeu et de la valeur absolue 
du coefficient d'ondelettes quantify voisin nord (ou ouest). Dans ce cas, le 



module de codage entropique est preferentiellement agence pour trailer 
quatre situations distinctes. En I'absence de coefficient voisin nord, il 
identifie les valeurs de prediction d'esperance et de largeur de la fonction 
laplacienne aux valeurs de prediction d'esperance ouest et de largeur 
ouest. En I'absence de coefficient voisin ouest, il identifie les valeurs de 
prediction d'esperance et de largeur de la fonction laplacienne aux valeurs 
de prediction d'esperance nord et de largeur nord. En I'absence de 
coefficients voisins ouest et nord, il effectue un codage laplacien fixe. Enfin, 
en presence de coefficients voisins ouest et nord, il identifie les valeurs de 
prediction d'esperance et de largeur de la fonction laplacienne 
respectivement a la demi-somme des valeurs de prediction d'esperance 
nord et ouest et a la demi-somme des valeurs de prediction de largeur nord 
et ouest, 

son module d'estimation peut etre agence de maniere a estimer les 
parametres de prediction des premier et second jeux par regression. 

• Par exemple, le module d'estimation peut etre charge d'estimer les 
premiers parametres de prediction des premier et second jeux a^'aide 
d'une premiere regression obtenue en resolvant une premiere equation 
de type W[n+1,m] = e.W[n,m], oCi W[n,m] est la valeur d'un coefficient 
d'ondelettes et W[n+1,m] est la valeur du coefficient d'ondelettes .yoisin 
de W[n,m], et 6 est un premier parametre de prediction. Dans ce cas, le 
module d'estimation peut par exemple resoudre la premiere equation 
par la technique des moindres carres. Le premier parametre de 
prediction d'un premier ou second jeu etant alors egal au rapport entre 
les covariance et variance non normalises des valeurs des coefficients 
d'ondelettes associes au premier ou second jeu. 

• Par ailleurs, le module d'estimation peut etre charge, d'une part, 
d'estimer les deuxieme et troisieme parametres de prediction du premier 
jeu a I'aide d'une seconde regression obtenue en resolvant un premier 
systeme d'equations de type W[n,m] - e N .W[n,m-1] = a N + 
Pn W[n,m-1] , ou W[n,m] est la valeur d'un coefficient d'ondelettes, 
Wn, m .i est la valeur du coefficient d'ondelettes voisin nord de W n>m , et 9 N , 
a N et (B N sont respectivement des premier, deuxieme et troisieme 
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parametres de prediction du premier jeu, et d'autre part, d'estimer les 
deuxieme et troisieme parametres de prediction du second jeu a I'aide 
d'une seconde regression obtenue en resolvant un second systeme 
d'equations de type W[n.m]-e 0 .W[n-1,m] =a 0 + po W[n-1,m] , ou 
W n . 1>m est la valeur du coefficient d'ondelettes voisin ouest de W n . m , et 
80, do et p 0 sont respectivement des premier, deuxieme et troisieme 
parametres de prediction du second jeu. Dans ce cas, ie module 
d'estimation peut par exemple resoudre les premier et second systemes 
d'equation a I'aide de la technique des moindres carres, 

- son module de codage entropique peut etre agence, lorsqu'un coefficient 
d'ondelettes quantifie voisin presente une valeur nulle, de maniere a le 
substituer par le resultat d'un produit entre le pas de quantification, utilise 
par le module de quantification, et une constante choisie, strictement 
superieure a zero et inferieure ou egale a un (1), et preferentiellement 
egale a 1/3, 

- son module de codage entropique peut etre agence de maniere a proceder 
a une discretisation de la fonction laplacienne apres avoir effectue les 
predictions d'esperance et de largeur. En variante, le module de codage 
entropique peut §tre agence de maniere a proceder au codage par iteration 
apres avoir effectue les predictions d'esperance et de largeur. 

L'invention porte egalement sur un compresseur/decompresseur de 
donnees d'image numeriques, ou codec, comprenant un dispositif de 
compression du type de celui pr6sent§ ci-avant. 

D'autres caracteristiques et avantages de l'invention apparaTtront a 
I'examen de la description detaillee ci-apres, et des dessins annexes, sur 
lesqueis : 

- la figure 1 illustre schematiquement, sous la forme d'un diagramme bloc 
fonctionnel, un dispositif de compression de donnees d'image numeriques 
selon l'invention, et 

- la figure 2 illustre une decomposition d'image en sous-bandes de niveaux 
de resolution differents. 

Les dessins annexes pourront non seulement servir a completer 
l'invention, mais aussi contribuer a sa definition, le cas echeant. 



L'invention a pour objet de permettre la compression et la 
decompression de donnees d'images numeriques par transformee en 
ondelettes. 

La compression (ou le codage) par transformee en ondelettes 
consiste classiquement a decomposer une image I en sous-bandes W jik dites 
LH, HL, HH et LL, puis a transformer les donnees d'image numeriques l[n,m], 
qui definissent I'intensite et la couteur des pixels de I'image l f en ondelettes, et 
plus precisement en un tableau [n,m] de coefficients d'ondelettes, puis a 
quantifier les coefficients d'ondelettes avant de les coder de facon entropique. 

n et m sont ici des indices entiers qui varient dans les limites de 
I'image : n designe la colonne et appartient a 1'intervalle [0, n max - 1] et m 
designe la ligne et appartient a I'intervalle [0, m max - 1]. Une sous-bande W j>k 
est done definie par un tableau bidimensionnel de coefficients d'ondelettes 
Wj, k [ri,m] oD les indices n et m varient dans des sous-intervalles [0, n,^ / 2 j - 1] 
et [0, m m ax / 2 j - 1]. On note alors N = n max . m max le nombre de pixels de 
I'image I, et % = N / 2 2j le nombre de coefficients d'ondelettes dans une- 
sous-bande Wj, k . 

Les sous-bandes sont codees separement les unes des autres, la 
sous-bande passe bas etant codee avec un codage differentiel de type DPCM ' 
et les differences resultantes pouvant etre codees avec une methode similaire : 
a celle utilisee pour une sous-bande passe-haut. 

Par ailleurs, on designe par Wj. k [n,m-1] et W j>k [n-1 ,m] les voisins dits 
« nord » et « ouest » d'un coefficient d'ondelettes W J(k [n,m]. Les coefficients 
d'ondelettes a coder sont parcourus dans un ordre particulier (ordre 
lexicographique croissant sur le couple (n.m)) afin de garantir qu'un coefficient 
Wj, k [n,m] sera toujours cod§ apres ses voisins nord et ouest qui servent a la 
determination de sa valeur. 

La compression selon l'invention differe d'une compression par 
transformee en ondelettes classique notamment par le fait qu'elle comprend 
entre son etape de transformation (classique) et son etape de quantification 
(classique), une etape dans laquelle on estime, pour chaque coefficient 
Wj, k [n,m] de chaque sous-bande W j>k d'une image I, de premier et second 
jeux de parametres de prediction, associes respectivement a des directions 
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nord-sud et ouest-est, en fonction des valeurs des coefficients d'ondelettes de 
ses voisins nord W j>k [n,m-1] ou ouest W i(k [n-1,m]. 

Ces premier et second jeux de parametres de prediction sont ensuite 
utilises lors de I'etape de codage entropique pour determiner pour chaque 
coefficient d'ondelettes W j|k [n,m] de chaque sous-bande Wj, k de I'image I des 
valeurs de prediction de I'esperance (ou moyenne) p et de la largeur a d'uhe 
fonction (ou distribution) de type laplacienne representative de sa densite de 
probability en fonction d'un coefficient d'ondelettes quantifie voisin (nord ou 
ouest) : 



ou c represente un coefficient d'ondelettes Wj |k [n,m], c' represente soit le 
coefficient d'ondelettes voisin nord Wj, k [n,m-1] de c, soit le coefficient 
d'ondelettes voisin ouest W j(k [n-1,m] de c, et u(c') et o(c') represented 
respectivement I'esperance et la largeur en fonction du coefficient 
d'ondelettes voisin c' considere. 

II est ici rappele que I'esperance p et la largeur a d'une distribution 
laplacienne de densite P sont definies par les relations suivantes : 

fi = jxP(x)dx 

cr = j\x-/x\P(x)dx 

On va maintenant decrire plus en detail les differentes etapes de la 
compression selon I'invention. 

Comme indique ci-avant, la compression commence par une etape 
de transformation des pixels de I'image I en coefficients d'ondelettes 
Wj,k[n,m]. Cette etape de transformation etant tout a fait classique, elle ne 
sera pas decrite ici. Pour la mettre en oeuvre on peut par exemple utiliser la 
transformee en ondelettes bi-orthogonales 7-9 de Daubechies/Antonini. Cette 
transformee n'est pas rigoureusement orthogonale, mais elle I'est quasiment, 
ce qui permet de conserver Penergie des images a quelques pour cent pres. 

A la fin de cette etape de transformation, on dispose done pour une 
sous-bande donnee W ]ik de I'ensemble de ses coefficients d'ondelettes 



p(c c') = 




2o-(c') 
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W j>k [n,m], 

L'etape de transformation est suivie d'une etape d'estimation de 
parametres de prediction pour chaque coefficient d'ondelettes W jik [n,m] de 
chaque sous-bande W jik . Dans ce qui suit, et afin de simplifier les notations on 
ne precisera plus les indices j et k des sous-bandes. Par consequent, on 
designera par W[n,m] un coefficient d'ondelettes appartenant a une sous- 
bande W jik [n,m], qu'il s'agisse d'une sous-bande LH, HH, ou HL, ou encore LL 
(par exemple codee en DPCM par un parcours en boustrophedon). 

Comme indique precedemment, l'etape d'estimation est destinee a 
determiner des premier et second jeux de parametres de prediction qui 
permettront ulterieurement (lors de I'etape de codage entropique) de predire 
les valeurs de I'esperance (ou moyenne) u et de la largeur a d'une fonction 
(ou distribution) laplacienne representative de la densite de probability d'un 
coefficient d'ondelettes quantifie compte tenu de la valeur prise par le 
coefficient d'ondelettes quantifie de I'un au moins de ses voisins nord et 
ouest. 

L'esperance u et la largeur a sont definies par rapport aux 
parametres de prediction des premier et second jeux. Par exemple, chaque: 
jeu comporte trois parametres de prediction 6 (premier), a (deuxieme) et @: 
(troisieme), et les relations liant ces parametres aux esperance u et largeur or. 
prennent les formes suivantes : 

M = G.c' (1) 

a = a + (3.c' (2) 

oD, comme indique precedemment c' designe le contexte representant soit le 
coefficient d'ondelettes voisin nord W[n,m-1] de c (ou W[n,m]), soit le 
coefficient d'ondelettes voisin ouest W[n-1 ,m] de c. 

En fait, comme on le verra plus loin, on utilise deux predictions 
d'esperance et de largeur selon le voisin considere (nord ou est) : une 
prediction nord porte sur I'esperance nord u N et la largeur nord o N , et une 
prediction ouest porte sur I'esperance ouest u N et la largeur ouest o N . Le 
premier jeu est alors associe a la prediction nord et comporte les parametres 
9n, a N et p N , tandis que le second jeu est associe a la prediction ouest et 
comporte les parametres 9 0 , a Q et p 0 - 
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Les relations (1) et (2) presentees ci-dessus, a titre d'exemple, 
definissent des regies de caicui qui permettent d'estimer simplement ies 
parametres G, a et p (nord et ouest) avec des regressions. 

Preferentiellement, on commence par estimer G N et 6 0 par une 
premiere regression. Celle-ci s'obtient en resolvant une equation du type : 

C n +1 = 6.C 

En utilisant, par exemple, ia technique des moindres carres, on 
obtient alors le premier parametre 8 pour chaque jeu (nord ou ouest) : 

£ = 

n 

La relation ci-dessus peut etre reecrite comme suit, pour les jeux nord 
et ouest : 

N Var 

e = Cov ° 

0 Var 

ou Cov N , Cov 0 et Var sont respectivement les covariances Nord et Ouest et la 
variance, non normalisees, donnees par les relations suivantes : 



Cov n = mty[n, m -l] 

n,m 



Puis, on procede a I'estimation des parametres de prediction nord a N 
et (3 N par une seconde regression. Celle-ci s'obtient en resolvant un premier 
systeme d'equations du type : 

W[n,m]-e N .W[n,m-1] =a N + (3 N W[n,m-1] (3) 

En imposant alors les definitions suivantes : 

^n= Zl (4) 



M l2 =M 2l =J]\ W h,™-l] (5) 
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Mn=Y,\w[n,m-l] 2 (6) 

ilytH 

Yi^Y^l n M-& N W[n t m-l\ (7) 



,M 21 M 2 JUJ IjJ 1 ; 



on obtient les parametres de prediction a N et p N en resolvant le systeme 
JUJ 

Les parametres de prediction ouest ao et po sont calcules de la 
m§me maniere, en resolvant un second systeme d'equations du type : 
W[n,m]-9 0 .W[n-1,m3 = a 0 + po W[n-1,m] (10) 

Cela revient a resoudre le systeme defini par la relation (9), en 
utilisant des coefficients definis par les relations (4) a (8) dans lesquelles e N 
est remplace par Go et W[n,m-1] est remplace par W[n-1 ,m], et en remplagant 
a N par a 0 et par po- 

Une fois que Ton a determine les parametres de prediction nord (0n, 
ct N et p N ) et ouest (0 O , ao et Po) de chaque coefficient d'ondelettes W Jtk [n,m] 
d'une sous-bande Wj, k , on passe a I'etape de quantification. 

Cette etape de quantification etant tout a fait classique, elle ne sera 
pas decrite ici. Pour la mettre en ceuvre on peut par exemple utiliser un 

quantificateur scalaire standard du type Q(W[n,m])= w ^ m ^ ou y es t j a 

largeur du pas de quantification. 

II est simplement rappele que I'etape de quantification consiste a 
quantifier les coefficients d'ondelettes Wj, k [n,m] de chaque sous-bande W jtk a 
partir d'un pas de quantification T soit fixe, soit determine par un module de 
planification. 

En raison de la nature meme des precedes de compression, le d6bit 
de sortie (du dispositif de compression) est directement proportionnel a 
I'entropie du signal d'entree (donnees d'images numeriques a compresser). 
Or, la quantit6 d'information a transmettre varie generalement au cours du 
temps puisque les caracteristiques du signal d'entree varient elles-memes 
avec le temps. Par consequent, dans le but d'optimiser I'utilisation du canal de 
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transmission, on minimise les fluctuations de debit en sortie du dispositif de 
compression en effectuant une regulation de debit. Cela est meme 
indispensable lorsque le canal de transmission necessite un debit constant. 

Pour permettre une telle regulation de debit en temps reel, il est 
imperatif de mettre en ceuvre des algorithmes de planification de debit rapides 
(exempts d'optimisation sous contraintes et exempts d'iterations) de maniere 
a ne pas ralentir I'acquisition des images lors du defilement du satellite ou, 
dans le cadre du codage video, pour ne pas generer d'effet de saccade. 

Deux types de planification de debit ont ete proposes a cet effet : la 
planification avec allocation de debit et la planification de debit pure. 

La planification avec allocation de debit consiste a laisser le module 
de planification, d'une part, choisir independamment les pas de quantification 
des differentes sous-bandes de sorte que le debit total soit egal a une 
consigne de debit, et d'autre part, choisir la facon dont ce debit doit etre 
reparti dans les differentes sous-bandes. 

La planification de debit pure consiste a lier les pas de quantification 
des differentes sous-bandes (en general ils sont identiques). Dans ce cas, le 
module de planification se contente de choisir le pas de quantification unique 
pourtoutes les sous-bandes qui permettra d'atteindre une consigne de debit 
fixee. 

Par ailleurs, afin de rel§cher la contrainte de debit fixe on peut 
controler la qualite de la compression (et done de la decompression) pour 
imposer une Constance de la qualite des images en sortie du dispositif de 
compression au gre du defilement du satellite ou de 1'arrivee des images en 
codage video. 

A cet effet, il a ete propose un procede de planification de qualite 
consistant a determiner un pas de quantification uniforme T permettant 
d'obtenir une qualite de codage donnee (generalement fixee en dB). Ce 
procede consiste tout d'abord a calculer la distorsion totale D correspondant a 
un pas de quantification donne T, puis a determiner par dichotomie le pas de 
quantification T permettant d'atteindre une consigne de distorsion donnee D. 

D'une facon plus precise, dans le domaine de I'allocation de debit ou 
du controle de qualite il existe principalement trois types de procede. 
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Un premier type concerne les precedes qui realisent une optimisation 
sous contraintes par une minimisation d'un lagrangien sous la forme J = D+ 
AR (ou D est la distorsion et R ie debit), avec pour contrainte R = Ro (ou Ro 
est un debit de consigne). 

Un deuxieme type concerne les precedes fondes sur une 
modelisation de la chaTne de compression pour etablir une regie qui lie le 
debit R au pas de quantification global T ou bien la distorsion D au debit R. 

Ces deux precedes ne peuvent cependant pas §tre utilises lorsque 
Ton a besoin de travailler avec des decompositions non uniforme, comme par 
exemple lors d'une re-decomposition de sous bandes pour la prise en compte 
de la FTM (Fonction de Transfert de Modulation) en vue d'une deconvolution 
ou d'une allocation par sous-bandes dans le cas d'images multicomposantes. 

Un troisieme type concerne ies precedes fondes sur une modelisation 
de la chaTne de compression pour etablir une regie qui lie le debit R au pas de 
quantification local T a chacune des sous-bandes ou bien la distorsion D au 
debit R. Un precede de ce type, adapte a I'allocation de debit par sous- 
bandes a I'aide d'un codeur arithmetique d'ordre 1 , est notamment decrit dans, 
le document FR 02/05724. i 

Un autre but de I'invention est d'ameliorer le resultat procure par lei 
compresseur (ou codeur) decrit dans le document precite en etablissant une 
regie optimale de pianification du debit ou de la qualite d'image d'un 
compresseur par transformee en ondelettes gr§ce au gain de codage 
engendre par la prise en compte des dependances entre coefficients 
d'ondelettes intra-bande. 

Pour atteindre ce but, il est tout d'abord necessaire de modeliser le 
comportement d'une sous-bande en termes de debit R et de distorsion D, en 
fonction des parametres de prediction estimes 9, a et p de la fonction 
laplacienne conditionnelle, determines lors de Petape d'estimation presentee 
ci-avant. 

Pour effectuer cette modelisation, il est possible d'utiliser 
simultanement, pour chaque sous-bande, les deux directions Nord-Sud et 
Ouest-Est, ou en d'autres termes les premier (9 N , a N , 3n) et second (9 0 , a 0 , 
Po) jeux de parametres de prediction. Mais, cela imposerait d'effectuer la 
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carte d'une fonction de grande dimension, et de plus, on peut montrer que le 
codage d'un coefficient d'ondelettes connaissant un seul coefficient 
d'ondelettes non quantifie voisin (nord ou ouest) donne un resultat de codage 
semblable a un codage utilisant les coefficients d'ondelettes non quantifies 
des voisins nord et ouest. Par consequent, dans ce qui suit on se place dans 
une situation dans laquelle on n'utilise qu'un seul coefficient d'ondelettes non 
quantifie voisin (nord ou ouest). ' 

Dans cette situation « monodirectionnelle », on peut definir une 
chame de dependance, dite « de Markov », par la relation suivante : 

P(X n Xn-i) = L _ e H^-^-,|/(«+^|^-,|) 

On peut alors montrer que si pour un jeu de parametres de prediction 
donnes 9, a et p, I'esperance et la variance de la variable aleatoire X„ restent 
bornees quand n croTt, alors la suite de variables aleatoires est ergodique. 

Pour calculer le debit R et la distorsion D d'une telle suite, pour un 
pas de quantification T donne, on peut utiliser une methode de Monte-Carlo 
consistent a faire des tirages d'une sequence de la chaTne de Markov, puis a 
estimer les esperances de debit R et de distorsion D par symbole en 
effectuant une moyenne sur des termes consecutifs de la sequence. 

Ces calculs de debit R et de distorsion D presentent un sens lorsque 
la chaTne de Markov est ergodique, et notamment des que sa largeur (ou 
variance) est bornee. 

Les tirages de la chaTne de Markov peuvent etre faits de la maniere 
suivante. On calcule tout d'abord par recurrence Xn +1 en fonction de X n par 
une formule du type : 

X n+ i = S n (a + p Xn )logU n + 9X n 

ou S n et U n sont des variables aleatoires independantes, chaque S n 
vaiant +1 et -1 avec des probabilites respectives 1/2 et 1/2, et U n ayant une 
densite uniforme sur rintervalie [0,1]. 

On doit ensuite estimer des courbes de debit et de distorsion. Pour ce 
faire, on doit tout d'abord estimer les relations qui lient le pas de quantification 
T au debit R et a la distorsion D, pour une gamme de parametres de 
prediction 9, a et p pertinente. 
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Par definition, la gamme des valeurs pertinentes est celle pour 
iaqueiie ia chatne de Markov est ergodique. On peut remarquer que 
I'ergodicite ne depend pas du parametre de prediction a qui est un parametre 
d'echelle. En effet, si une chaine de variables aleatoires X n est une chatne de 
Markov pour les parametres de prediction 9, a et p, alors la chatne KK n est 
une chatne de Markov pour les parametres 9, Act et p. On peut done se limiter 
au cas a = 1 . 

Par ailleurs, il est naturel de limiter les valeurs de (3 a I'intervalle [0,1], 
puisque p ne peut pas etre negatif. Si p est plus grand que 1, il est 
vraisemblable que la variance de la distribution associee tende vers I'infini 
(par un etalement progressif). De meme, 9 doit presenter une valeur absolue 
inferieure a 1 , sans quoi le processus risquerait d'avoir un comportement de 
marche aleatoire, lequel n'est pas ergodique. En revanche, le domaine des 
parametres de prediction 9 licites (les autres parametres de prediction etant 
fixes) est necessairement symetrique, puisque si une suite Xn est de 
parametres 9, a et p, la suite (-1)% est de parametres -9, a et p. Par un, 
raisonnement analogue, on peut meme montrer que le domaine en question „ 
est un intervalle de la forme ]-9o, 9o[. 

Faute de pouvoir definir clairement par des methodes analytiques un ;,;, 
domaine des couples (P, 9) pour lequel la chatne de Markov est ergodique, on. , 
en definit un sous-ensemble pour lequel la largeur estimee de la variable 
aleatoire E( X„ ) est bornee par une limite fixee. 

Lorsque Ton trace la carle de la variable aleatoire E(Xn) en fonction 
des parametres de prediction p et 9, on constate que I'estimation d'une 
esperance devient de plus en plus coGteuse (en nombre d'iterations) au fur et 
a mesure que Ton s'approche de la non-ergodicite. Par ailleurs, en tracant des 
sections de cette fonction de deux variables et de son inverse, on peut 
cpnstater que son inverse est sensiblement lineaire en p : 

E(X„)=: 

ou A et B sont des fonctions a valeurs positives. 
Par consequent, pour chaque section de cette fonction a deux 
variables, a 9 constant, la pente et I'ordonnee a I'origine peuvent etre 
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estimees par une Agression tres stable. Par exemple, les fonctions A(8) et 
B(9) peuvent etre parametrees par decomposition dans une base de 
Deslauriers-Dubuc. II est done globalement possible d'estimer les limites 
d'ergodicite avec une grande precision. 

Le parametrage des fonctions A(9) et B(9) au voisinage de un (1) sont 
plus precis si on decompose dans une base de Deslauriers-Dubuc non pas 
les fonctions A(9) et B(8), mais plutdt A(9) 2 et B(9) 2 . En effet les graphes de 
ces fonctions semblent montrer qu'elles ont un zero de la forme (1 - 8) 1/2 en 
un (1) pour A, et pour une valeur 9 > 1 pour B. La decomposition dans la base 
de Deslauriers-Dubuc de chaque fonction coQte neuf coefficients. 

Si I'on appelle R(8, a, p, T) la fonction de debit d'une sous-bande en 
fonction des parametres de prediction 8, a et p et du pas de quantification T, 
on peut la definir par la relation suivante : 

R(9, a, 3, T) = E(-P(Q(X n ) Xn-,) log 2 P(Q(Xn) VO) 

ou Xn est une chaine de Markov laplacienne de parametres 9, a et p, 
et Q un quantificateur de pas T. 

Le debit R peut etre estime par des moyennes de Monte-Carlo 
calculees sur des tirages de la chaine de Markov, done par moyennage de la 
suite : 

Rn = -P n log 2 P n 
oil 

p n =~ few = Q{x n ) — i_ e -i-^- 1 |/(^^,i>^ 

Ces moyennes peuvent §tre estimees pour des parametres de 
prediction 9 et p qui sont dans le domaine d'ergodicite identifie plus haut. On 
peut utiliser I'homogeneite de la fonction de debit R pour reduire la 
dimensionnalite du probleme, du fait que R(8, a, p, T) = R(9, 1, p, T/ a). 

Si I'on appelle maintenant D(9, a, p, T) la fonction de distorsion de la 
chatne de Markov, qui servira a modeliser la distorsion D d'une sous-bande 
quantifiee en fonction des parametres de prediction 8, a et p et du pas de 
quantification T, on peut la definir par la relation suivante : 

D(9,a,p,T)=E( Q(X n )-X n 2 ) 

ou X n est de nouveau la chaine de Markov laplacienne de parametres 
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9, a et 0, et Q un quantificateur de pas T. 

Ici encore, la distorsion D peut etre estimee par des moyennes de 
Monte-Carlo de la suite D n = X„ - Q(X n ) . II faut en effet connattre la 
distribution stationnaire associee a la chaTne de Markov. Or, il n'en existe pas 
d'expression explicite, de sorte qu'il faut recourir a une methode de Monte- 
Carlo. 

Ici encore, I'homogeneite permet de reduire la dimension du 
probteme : 

D(6, a,p,T) = a 2 D(9, 1,p,T/a) 

Etant donne que le domaine d'ergodicite contient seulement trois 
points singuliers (ou coins), on peut utiliser un maillage triangulaire du 
domaine d'ergodicite concerne et utiliser des fonctions splines triangulaires 
pour representer les fonctions de distorsion D et de debit R sur ce domaine 
d'ergodicite. Cependant, on peut montrer que le coin superieur du domaine 
est un point ou les fonctions D et R peuvent avoir de grandes variations, et 
qu'il est preferable d'utiliser un maillage plus fin autour de ce point, comme 
par exemple un maillage carre projete dans le domaine d'ergodicite. 

Pour 9, p dans le domaine d'ergodicite maille, a = 1 et la variable T> 
(pas de quantification) variant entre et 1000, des estimations de la 
distorsion D et du debit R peuvent §tre faites par la methode de Monte-Carlo 
sur des grandes series (typiquement quelques millions de tirages). On obtient 
ainsi une carte de D et R sur un domaine d'ergodicite qui couvre une gamme 
de parametres tres large. 

Ensuite, la fonction inverse (qui associe un pas de quantification T a 
un debit R) peut etre calculee et decomposee dans une base de Deslauriers- 
Dubuc, par inversion de I'application 9, p, T R en 9, p, R T. 

Dans cette inversion, il faut prendre un parametrage en R soit 
lineaire, soit logarithmique, afin de representer fidelement la fonction T(R) 
dans des regions oD eile varie tres brutalement. 

On peut ensuite effectuer la prediction des relations entre le debit R, 
le pas de quantification T et la distorsion D pour une sous-bande donnee. 

Comme on I'a vu prec§demment, pour chaque sous-bande la 
prediction d'un coefficient d'ondelettes se fait preferentiellement a I'aide de 
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deux coefficients d'ondelettes voisins (nord et ouest), alors que les modeles 
de debit R et de distorsion D presentes ci-avant n'utilisent qu'un contexte a un 
seul coefficient d'ondelettes voisin. 

Pour associer un modele monodimensionnei a une sous-bande, il 
faut alors reduire la paire de premier et second jeux de parametres de 
prediction, associes a deux fonctions laplaciennes marginales, a un unique 
jeu. La reduction utilisee pour la planification de debit est par exemple la 
suivante : 

a = ^a N a Q 
P 2VT" 
2 

On peut voir cette reduction comme un moyennage soit arithm&ique, 
soit geometrique. Seul 0 est en plus encore contracts d'un facteur 2, cense 
modeliserla reduction d'incertitude liee a I'usage d'un contexte double. 

Le choix du facteur 2 peut atre justifie par une approche bayesienne 
simple : deux mesures, bruitees par des bruits independants combines, 
peuvent §tre combinees (par moyennage) en reduisant I'ecart-type du bruit de 
mesure d'un facteur 2. Le choix de cette convention est certes arbitraire, 
mais il permet d'obtenir des resultats raisonnables de prediction pour des 
images classiques. Bien entendu, d'autres choix peuvent §tre envisages pour 
affiner cette reduction. 

Compte tenu de ce qui a ete precise ci-avant, il est desormais 
possible d'elaborer simplement des algorithmes de planification. 

Ainsi, dans le cas d'une planification de debit pure, consistant a 
determiner un pas de quantification uniforme T permettant d'atteindre un debit 
de consigne global pour une image, il suffit de calculer pour chaque sous- 
bande une prediction de debit en fonction d'un pas de quantification T, soit 
done un debit total R. Pour ce faire, on determine par dichotomie la valeur du 
pas de quantification T qui rend le debit R, correspondant a T, egal a la 
consigne de debit. 

Dans le cas d'une planification par allocation de debit, on cherche a 
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determiner librement un pas de quantification T pour chaque sous-bande, afin 
de minimiser la distorsion D pour une consigne de debit globale R. Pour ce 
faire, on decoupe le debit total R en trongons de petite taille, et on distribue 
ces trongons de fagon iterative entre les differentes sous-bande, en donnant a 
chaque fois un trongon a ia sous-bande qui reduit le plus la distorsion globale 
D. Comme indique precedemment, un exemple de mise en ceuvre de cette 
methode est decrit dans le document FR 02/05 724. 

Dans le cas d'une planification de qualite, on cherche a determiner un 
pas de quantification uniforme T permettant d'obtenir une qualite de codage 
donnee (fixee en dB). Pour ce faire, on calcuie la distorsion totale D qui 
correspond a un pas de quantification donne T, puis on determine par 
dichotomie le pas de quantification T qui permet d'atteindre la consigne de 
distorsion donn§e D. 

A la fin de I'etape de quantification, avec ou sans planification, on 
dispose pour une sous-bande donnee W itk de I'ensemble de ses coefficients 
d'ondelettes quantifies Q(W j(k [n,m]). 

On passe alors a I'etape de codage entropique dit « contextuel ». 
Cette etape se decompose en deux phases, I'une dediee a la prediction des 
esperances p et largeurs a, a partir, notamment, des premiers et seconds 
jeux de parametres de prediction determines lors de I'etape d'estimation,. et 
l-autre dediee au codage entropique proprement dit. 

Le mode de prediction de Pesperance p et de ia largeur a associees a 
un coefficient d'ondelettes quantifie depend de son contexte, ou en d'autres 
termes des coefficients d'ondelettes voisins (lorsqu'ils existent). Quatre 
situations peuvent etre rencontrees. 

En Tabsence de coefficient d'ondelettes voisin nord (W[n,m-1]), c'est- 
a-dire lorsque le coefficient d'ondelettes (W[n,1]) appartient a la premiere 
ligne de la sous-bande concernee, on identifie les valeurs de prediction 
d'esperance p et de largeur a de la fonction laplacienne (P) aux valeurs de 
prediction d'esperance ouest po et de largeur ouest a 0l soit p = po et a = a 0 . 

En I'absence de coefficient d'ondelettes voisin ouest (W[n-1,m]), 
c'est-a-dire lorsque le coefficient d'ondelettes (W[l,m]) appartient a la 
premiere colonne de la sous-bande concernee, on identifie les valeurs de 
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prediction d'esperance p et de largeur a de la fonction lapiacienne (P) aux 
valeurs de prediction d'esperance nord M n et de largeur nord o M , soit p = p N et 
a = a N . 

En I'absence de coefficients d'ondelettes voisins ouest (W[n-1,m]) et 
nord (W[n,m-1]), c'est-a-dire lorsque le coefficient d'ondelettes (W[1,1]) 
appartient a la premiere ligne et a la premiere colonne de la sous-bande 
concernee, on effectue un codage laplacien fixe, c'est-a-dire un codage 
laplacien d'ordre 0 fonde sur la prise en compte d'une density de probability 
non conditionnelle. 

Enfin, en presence de coefficients d'ondelettes voisins ouest et nord, 
on identifie les valeurs de prediction d'esperance p et de largeur a de la 
fonction lapiacienne (P) respectivement a la demi-somme des valeurs de 
prediction d'esperance nord p N et ouest Mo et a la demi-somme des valeurs 
de prediction de largeur nord o N et ouest cr 0 , soit p = [(p N + Mo ) / 2] et a = [(a N 
+ a 0 )/2]. 

Les valeurs de p N , p 0 , o N et o 0 sont determinees a partir des relations 
(1) et (2) qui sont mentionn6es de nouveau ci-apres pour les predictions nord 
et ouest : 

Pn = e N .Q(W[n,m-1]) (1a) 

Po = 9 0 .Q(W[n-1,m]) (1b) 

a N = a N + 0 N . Q(W[n,m-1]) (2a) 

a o = a o + p 0 . Q(W[n-1,m]) (2b) 
oD Q(W[n,m-1]) et Q(W[n-1,m]) sont les coefficients d'ondelettes quantifies 
des voisins nord et ouest d'un coefficient d'ondelettes quantifie Q(W[n,m]). 

II est important de noter que le calcul de la largeur o est 
generalement biaisee lorsque la valeur du coefficient d'ondelettes quantifie 
Q(W[n,m-1]) ou Q(W[n-1 ,m]) associe est nulle. 

Dans ce cas, en constatant que Q(W[n,m]) represente I'esperance 
E(W[n,m] Q(W[n,m])) de W[n,m] connaissant sa valeur quantifiee, on peut 
en deduire la relation suivante lorsque Q(W[n,m]) = 0 : 

E( W[n,m] Q(W[n,m])) = T/2 Q(W[n,m]) 
ou T est le pas de quantification choisi a I'etape de quantification par le 
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quantificateur (suppose ici etre du type dit « a double boTte 0 »). 

Et, partant de cette deduction, on peut alors remplacer dans les 
relations (1a), (1b), (2a) et (2b) Q(W[n,m-1]) et Q(W[n-1,m]), lorsqu'ils sont 
nuls, par une constante £T choisie. £ est une constante strictement superieure 
a zero et inferieure ou egale a 1. Une valeur de £ 6gale a 1 / 2 est en theorie 
optimale, mais, dans le cas d'une quantification a haute resolution et de 
fonctions laplaciennes tres pentues dans la bolte 0, il est preferable de choisir 
Regale a 1/3. 

Une fois que Ton est en possession de I'esperance u et de la largeur 
o associees a chaque coefficient d'ondelettes quantifie d'une sous-bande, on 
peut passer a la phase de codage entropique. 

Pour ce faire, il suffit de discretiser une laplacienne qui n'est pas 
necessairement centree dans la boTte 0 du quantificateur. La discretisation 
consiste ici a echantillonner la laplacienne continue. 

La decompression de donnees d'images numeriques, compressees 
au moyen du procede de compression qui vient d'etre decrit, consiste a 
enchaTner les quatre etapes dudit procede dans I'ordre inverse. Par 
consequent, le procede de decompression ne sera pas decrit ci-apres. 

L'invention porte egalement sur un dispositif de compression (ou 
codeur) de donnees d'images numeriques par transformee en ondelettes. Ce 
dispositif ayant pour objet de mettre en ceuvre le procede de compression 
decrit ci-avant, seuls ses modules fonctionnels seront presentes ci-apres en 
reference a Punique figure. 

Le dispositif de compression D, selon l'invention, comprend un 
module de traitement PM decomposable en au moins quatre modules. 

Un premier module M1 est alimente en donnees d'image numeriques 
l[n,m] definissant des images I et est charge, d'une part, de decomposer 
chaque image I sur une base d'ondelettes. Cette transformation en sous- 
bandes W jik (dites LH, HL, HH et LL comme indique sur la figure 2) et sur 
piusieurs niveaux de resolution j. Chaque sous-bande est representee par un 
tableau [n,m] de coefficients d'ondelettes W j|k [n,m]. II s'agit d'un module 
traditionnel, parfaitement connu de I'homme de I'art, mettant en ceuvre, par 
exemple, la transformee en ondelettes bi-orthogonales 7-9 de 
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Daubechies/Antonini. 

Un deuxieme module M2 est alimente en coefficients d'ondelettes 
Wj, k [n,m] par le premier module de transformation M1. II est charge d'estimer, 
pour chaque coefficient d'ondelettes W jik [n,m] de chaque sous-bande 
d'ondelettes W Jik d'une image I, des premier (6 N , a N et (3 N ) et second (9 0 , a Q et 
Po) jeux de parametres de prediction, associes respectivement a des 
directions nord-sud et ouest-est, en fonction des valeurs des coefficients 
d'ondelettes de ses voisins nord W j>k [n,m-1] ou ouest W jik [n-1,m]. Pour ce 
faire, le deuxieme module d'estimation M2 met en ceuvre par exemple la 
technique d'estimation a double regression decrite precedemment. 

Un troisieme module M3 est alimente en coefficients d'ondelettes 
Wj.ktn.m] par le premier module de transformation M1 et est charge de les 
quantifier pour delivrer des coefficients d'ondelettes quantifies Q(W j>k [n,m]). II 
s'agit d'un module de quantification scalaire (ou quantificateur) traditionnel, 
parfaitement connu de I'homme de Tart. 

Un quatrieme module M4 est alimente, d'une part, en coefficients 
d'ondelettes quantifies Q(Wj ik [ n> m]) par le troisieme module de quantification 
M3, et d'autre part, en premier et second jeux de parametres de prediction (6, 
a et 0) par le deuxieme module d'estimation M2. II est charge de coder de 
facon entropique et contextuelle les coefficients d'ondelettes quantifies 
Q(W j>k [n,m]) a partir des valeurs des parametres de prediction nord (G Nl a N et 
|3 N ) et ouest (9 0 , a 0 et 0 O ) des premier et second jeux. 

Pius precisement, le quatrieme module d'estimation M4 est agence 
de maniere a determiner pour chaque coefficient d'ondelettes quantifies 
Q(W j , k [n 1 m]), les valeurs de prediction d'esperance p et de largeur a de la 
fonction laplacienne P associee, compte tenu des valeurs des coefficients 
d'ondelettes quantifies de ses voisins nord W j>k [n,m-1] et/ou ouest Wj, k [n-1,m] 
et des parametres de prediction nord (9 N , a N et f3 N ) et ouest (9 0 , a 0 et p 0 ) 
associes, puis a coder de facon entropique chaque coefficient d'ondelettes 
quantifie Q(W jik [n,m]) compte tenu des valeurs de prediction d'esperance u et 
de largeur a qu'il a determines. Pour ce faire, le quatrieme module de 
codage entropique M4 met en ceuvre par exemple les techniques de 
determination de prediction et de codage entropique contextuel decrites 
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prec6demment. 

Une fois le codage entropique termini, !e quatrieme module de 
codage entropique M4 delivre sur sa sortie des donnees d'images 
numeriques compressees DC. 

Le dispositif de compression D peut egalement comporter un 
cinquieme module M5 charge de la planification de debit et/ou de I'allocation 
de debit par sous-bandes et/ou de la planification de qualite, 
preferentiellement selon le procede de planification presente ci-avant. Ce 
cinquieme module M5 est couple au deuxieme module d'estimation M2 qui 
I'alimente en premier (0 N , ct N et p N ) et second (9 0 , a 0 et 0 O ) jeux de 
parametres de prediction et au troisieme module de quantification M3. 

L'invention porte egalement sur un dispositif de decompression (ou 
decodeur) mettant en ceuvre le procede de decompression precite. Ce 
dispositif de decompression ne sera pas decrit ci-apres dans la mesure ou il 
est sensiblement constitue des memes modules fonctionnels que le dispositif 
de compression D, mais avec un fonctionnement inverse (le module de 
codage etant ici alimente en donnees d'images numeriques compressees et 
le module de transformation delivrant en sortie des donnees d'images 
numeriques decompressees definissant I'image initiale). 

L'invention porte en outre sur un compresseur/decompresseur, ou 
codeur/decodeur, ou encore codec, comprenant un dispositif de compression 
et un dispositif de decompression selon l'invention. 

L'invention offre un procede de codage particulierement efficace 
assurant un gain de codage notable grace a I'exploitation de I'information 
mutuelle conjointement a la modelisation laplacienne conditionnelle. 

Par ailleurs, grace aux techniques d'estimation des parametres de 
prediction et de determination des predictions, particulierement simples en 
termes de calcul, les structures des codeur et decodeur sont simplifies. 

En outre, du fait que Ton dispose d'un quantificateur different pour 
chacune des sous-bandes, il est done possible de proceder a une allocation 
locale. 

Enfin, du fait que la regie d'allocation prend explicitement en compte 
les caracteristiques du dispositif de codage entropique d'ordre 1 selon 
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I'invention, ledit dispositif de compression est beaucoup plus efficace que 
ceux de I'art anterieur, notamment pour les applications en temps reel a haut 
debit d'entree (grandes images), en particulier dans le domaine spatial. Cette 
efficacite est encore renforcee du fait du faible nombre d'operations requises 
par le codage. 

L'invention ne se limite pas aux modes de realisation de precede de 
compression, precede de decompression, dispositif de compression, dispositif 
de decompression et codec decrits ci-avant, seulement a titre d'exemple, 
mais elle englobe toutes les variantes que pourra envisager I'homme de i'art 
dans le cadre des revendications ci-apres. 



27 



REVENDICATIONS 

1. Procede de compression de donnees d'images numeriques par 
transformee en ondelettes, comprenant une etape de transformation des 
donnees d'image numeriques en coefficients d'ondelettes repartis dans des 
sous-bandes, une etape de quantification desdits coefficients d'ondelettes, et 
une etape de codage entropique desdits coefficients d'ondelettes quantifies, 
caracterise en ce qu'il comprend, entre les etapes de transformation et de 
quantification, une etape d'estimation, pour chaque coefficient de chaque 
sous-bande de I'image, de premier et second jeux de parametres de 
prediction, associes respectivement a des directions dites « nord-sud » et 
« ouest-est », en fonction des valeurs des coefficients d'ondelettes de ses 
voisins nord ou ouest, et en ce qu'a I'etape de codage entropique, d'une part, 
on determine pour chaque coefficient de chaque sous-bande de I'image une 
prediction de I'esperance et de la largeur d'une fonction de type laplacienne, 
representative de sa densite de probability, en fonction desdits premier et/ou 
second jeux de parametres de prediction et du coefficient d'ondelettes 
quantifie de son voisin nord ou de son voisin ouest, et d'autre part, on code de 
facon entropique lesdits coefficients d'ondelettes quantifies a I'aide desdites 
esperances et largeurs associees determinees. 

2. Procede selon la revendication 1 , caracterise en ce qu'a I'etape de 
codage entropique on determine pour la fonction laplacienne de chaque 
coefficient d'ondelettes quantifie, lorsque cela est possible, des predictions 
d'esperance nord et ouest et des predictions de largeur nord et ouest, ladite 
prediction d'esperance nord, respectivement ouest, etant definie par le produit 
entre un premier parametre de prediction du premier, respectivement second, 
jeu et la valeur dudit coefficient d'ondelettes quantifie voisin nord, 
respectivement ouest, et ladite prediction de largeur nord, respectivement 
ouest, etant definie par la somme d'un deuxieme parametre de prediction du 
premier, respectivement second, jeu et du produit d'un troisieme parametre 
de prediction du premier, respectivement second, jeu et de la valeur absolue 
dudit coefficient d'ondelettes quantifie voisin nord, respectivement ouest 

3. Procede selon la revendication 2, caracterise en ce qu'a I'etape de 
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codage entropique, d'une premiere part, en I'absence de coefficient voisin 
nord on identifie les predictions d'esperance et de largeur de la fonction 
laplacienne aux predictions d'esperance ouest et de largeur ouest, d'une 
deuxieme part, en I'absence de coefficient voisin ouest on identifie les 
predictions d'esperance et de largeur de la fonction lapiacienne aux 
predictions d'esperance nord et de largeur nord, d'une troisieme part, en 
I'absence de coefficients voisins ouest et nord on effectue un codage 
laplacien fixe, et d'une quatrieme part, en presence de coefficients voisins 
ouest et nord on identifie les predictions d'esperance et de largeur de la 
fonction laplacienne respectivement a la demi-somme des predictions 
d'esperance nord et ouest et a la demi-somme des predictions de largeur 
nord et ouest. 

4. Procede selon I'une des revendications 1 a 3, caracterise en ce que 
Ton estime les parametres de prediction desdits premier et second jeux par 
regression. 

5. Procede selon la combinaison des revendications 2 et 4, caracterise 
en ce que I'on estime lesdits premiers parametres de prediction des premier 
et second jeux par une premiere regression obtenue par resolution d'une 
premiere equation W[n+1,m] = e.Wfn.m] oD W[n,m] est la valeur d'un 
coefficient d'ondelettes et W[n+1,m] est la valeur du coefficient d'ondelettes 
voisin de W[n,m], et 9 est un premier parametre de prediction. 

6. Procede selon la revendication 5, caracterise en ce que I'on resout 
ladite premiere equation par la technique des moindres carres, ledit premier 
parametre de prediction d'un premier ou second jeu etant egal au rapport 
entre les covariance et variance non normalisees des valeurs des coefficients 
d'ondelettes, associes au premier ou second jeu. 

7. Procede selon la revendications 2 en combinaison avec Tune des 
revendications 4 a 6, caracterise en ce que I'on estime lesdits deuxieme et 
troisieme parametres de prediction du premier jeu par une seconde 
regression obtenue par resolution d'un premier systeme d'equations W[n,m] 
- e N .W[n,m-1] = a N + (3 N W[n,m-1] , ou W[n,m] est la valeur d'un 
coefficient d'ondelettes, W n . m .i est la valeur du coefficient d'ondelettes voisin 
nord de W n , m , et 0 N , a N et p N sont respectivement des premier, deuxieme et 
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troisieme parametres de prediction du premier jeu, et en ce que Ton estime 
lesdits deuxieme et troisieme parametres de prediction du second jeu par une 
seconde regression obtenue par resolution d'un second systeme d'equations 
W[n,m] - 9 0 .W[n-1,m] = a 0 + po W[n-1,m] , ou W n . 1>m est !a valeur du 
coefficient d'ondelettes voisin ouest de W n , m , et 0 Ol a D et (3 0 sont 
respectivement des premier, deuxieme et troisieme parametres de prediction 
du second jeu. 

8. Procede selon la revendication 7, caracterise en ce que Ton resout 
lesdits premier et second systemes d'equation par la technique des moindres 
carres. 

9. Procede selon Tune des revendications 1 a 8, caracterise en ce 
qu'en presence d'un coefficient d'ondelettes quantifie voisin de valeur nulie, 
on le substitue par un produit entre le pas de quantification, utilise lors de 
Tetape de quantification, et une constante choisie, strictement superieure a 
zero et inferieure ou egale a un (1). 

10. Proced6 selon la revendication 9, caracterise en ce que iadite 
constante est choisie 6gale a 1/3. 

11. Procede selon Tune des revendications 1 a 10, caracterise en ce 
qu'a Tetape de codage entropique, apres avoir effectue les predictions 
desdites esperance et largeur, on procede a la discretisation de la fopction 
laplacienne. 

12. Procede selon Tune des revendications 1 a 10, caracterise en ce 
qu'a retape de codage entropique, apres avoir effectue les predictions 
desdites esperance et largeur, on procede audit codage par iteration. 

13. Dispositif (D) de compression de donnees d'images numeriques par 
transformee en ondelettes, comprenant des moyens de traitement (PM) 
munis, d'une premiere part, d'un module (M1) de transformation de donnees 
d'image numeriques en coefficients d'ondelettes repartis dans des sous- 
bandes, d'une deuxieme part, d'un module (M3) de quantification des 
coefficients d'ondelettes, et d'une troisieme part, d'un module (M4) de codage 
entropique des coefficients d'ondelettes quantifies, caracterise en ce que 
lesdits moyens de traitement (PM) comprennent en outre un module 
d'estimation (M2) agence pour estimer pour chaque coefficient de chaque 
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sous-bande de I'image, de premier et second jeux de parametres de 
prediction, associes respectivement a des directions dites « nord-sud » et 
« ouest-est », en fonction des valeurs des coefficients d'ondelettes de ses 
voisins nord ou ouest, et en ce que ledit module de codage entropique (M4) 
est agence, d'une part, pour determiner pour chaque coefficient de chaque 
sous-bande de I'image des valeurs de prediction de I'esperance et de la 
largeur d'une fonction de type laplacienne, representative de sa densite de 
probabilite, en fonction desdits premier et/ou second jeux de parametres de 
prediction et du coefficient d'ondelettes quantifie de son voisin nord ou de son 
voisin ouest, et d'autre part, de coder de facon entropique lesdits coefficients 
d'ondelettes quantifies a I'aide desdites esperances et largeurs associees 
determinees. 

14. Dispositif selon la revendication 13, caracterise en ce que ledit 
module de codage entropique (M4) est agence pour determiner pour la 
fonction laplacienne de chaque coefficient d'ondelettes quantifie, lorsque cela 
est possible, des valeurs de prediction d'esperance nord et ouest et des 
valeurs de prediction de largeur nord et ouest, ladite valeur de prediction 
d'esperance nord, respectivement ouest, etant definie par le produit entre un 
premier parametre de prediction du premier, respectivement second, jeu et la 
valeur dudit coefficient d'ondelettes quantifie voisin nord, respectivement 
ouest, et ladite valeur de prediction de largeur nord, respectivement ouest, 
etant definie par la somme d'un deuxieme parametre de prediction du 
premier, respectivement second, jeu et du produit d'un troisieme parametre 
de prediction du premier, respectivement second, jeu et de la valeur absolue 
dudit coefficient d'ondelettes quantifie voisin nord, respectivement ouest. 

15. Dispositif selon la revendication 14, caracterise en ce que ledit 
module de codage entropique (M4) est agence, d'une premiere part, en 
I'absence de coefficient voisin nord pour identifier les valeurs de prediction 
d'esperance et de largeur de la fonction laplacienne aux valeurs de prediction 
d'esperance ouest et de largeur ouest, d'une deuxieme part, en I'absence de 
coefficient voisin ouest pour identifier les valeurs de prediction d'esperance et 
de largeur de la fonction laplacienne aux valeurs de prediction d'esperance 
nord et de largeur nord, d'une troisieme part, en I'absence de coefficients 
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voisins ouest et nord pour effectuer un codage laplacien fixe, et d'une 
quatrieme part, en presence de coefficients voisins ouest et nord pour 
identifier les valeurs de prediction d'esperance et de largeur de la fonction 
laplacienne respectivement a la demi-somme des valeurs de prediction 

5 d'esperance nord et ouest et a la demi-somme des valeurs de prediction de 
largeur nord et ouest 

16. Dispositif selon Tune des revendications 13 a 15, caracteris§ en ce 
que ledit module d'estimation (M2) est agence pour estimer les parametres de 
prediction desdits premier et second jeux par regression. 

10 17. Dispositif selon la combinaison des revendications 14 et 16, 

caracterise en ce que ledit module d'estimation (M2) est agence pour estimer 
lesdits premiers parametres de prediction des premier et second jeux par une 
premiere regression obtenue par resolution d'une premiere equation 
W[n+1,m] = G.W[n,m] oCi W[n,m] est la valeur d'un coefficient d'ondelettes et 

is W[n+1,m] est la valeur du coefficient d'ondelettes voisin de W[n,m], et 9 est 
un premier parametre de prediction. 

*• • 

18. Dispositif selon la revendication 17, caracterise en ce quejedit 
module d'estimation (M2) est agence pour resoudre ladite premiere equation 
par la technique des moindres carres, ledit premier parametre de prediction 

20 d'un premier ou second jeu etant egal au rapport entre les covariance et 
variance non normalisees des valeurs des coefficients d'ondelettes, associes 
au premier ou second jeu. 

19. Dispositif selon la revendications 14 en combinaison avec Tune des 
revendications 16 a 18, caracterise en ce que ledit module d'estimation (M2) 

25 est agence, d'une part, pour estimer lesdits deuxieme et troisieme parametres 
de prediction du premier jeu par une seconde regression obtenue par 
resolution d'un premier systeme d'equations W[n,m] - e N .W[n,m-1] = a N + 
B N W[n,m-1] , oCi W[n,m] est la valeur d'un coefficient d'ondelettes, W n>m .i 
est la valeur du coefficient d'ondelettes voisin nord de W n>m , et 6 N , a N et 8 N 

30 sont respectivement des premier, deuxieme et troisieme parametres de 
prediction du premier jeu, et d'autre part, pour estimer lesdits deuxieme et 
troisieme parametres de prediction du second jeu par une seconde regression 
obtenue par resolution d'un second systeme d'equations W[n,m] - 0 o .W[n- 
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1.m] = a 0 + Po W[n-1,m] , oD W n . 1im est la valeur du coefficient 
d'ondeiettes voisin ouest de W n>m , et 6 0 , ct 0 et fr> sont respectivement des 
premier, deuxieme ettroisieme parametres de prediction du second jeu. 

20. Dispositif selon la revendication 19, caracterise en ce que ledit 
module d'estimation (M2) est agence pour resoudre lesdits premier et second 
systemes d'equation par la technique des moindres carres. 

21. Dispositif selon I'une des revendications 13 a 20, caracterise en ce 
que ledit module de codage entropique (M4) est agence, en presence d'un 
coefficient d'ondeiettes quantifie voisin de valeur nulle, pour le substituer par 
le resultat d'un produit entre le pas de quantification, utilise par le module de 
quantification, et une constante choisie, strictement superieure a zero et 
inferieure ou egale a 1 . 

22. Dispositif selon la revendication 21, caracterise en ce que ladite 
constante est choisie egale a 1/3. 

23. Dispositif selon I'une des revendications 13 a 22, caracterise en ce 
que ledit module de codage entropique (M4) est agence pour proceder a une 
discretisation de la fonction laplacienne apres avoir effectue les predictions 
desdites esperance et largeur. 

24. Dispositif seion I'une des revendications 13 a 22, caracterise en ce 
que ledit module de codage entropique (M4) est agence pour proceder audit 
codage par iteration apres avoir effectue les predictions desdites esperance et 
largeur. 

25. Codec de donnees d'image numeriques, caracterise en ce qu'ii 
comprend un dispositif de compression (D) selon I'une des revendications 13 
a 24. 
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